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ABSTRACT 

Molecular hydrogen emission is commonly observed in planetary nebulae. Images 
taken in infrared H2 emission lines show that at least part of the molecular emission is 
produced inside the ionised region. In the best-studied case, the Helix nebula, the H2 
emission is produced inside cometary knots (CKs), comet-shaped structures believed 
to be clumps of dense neutral gas embedded within the ionised gas. Most of the H2 
emission of the CKs seems to be produced in a thin layer between the ionised diffuse 
gas and the neutral material of the knot, in a mini photodissociation region (PDR). 
However, PDR models published so far cannot fully explain all the characteristics of 
the H2 emission of the CKs. In this work, we use the photoionisation code Aangaba 
to study the H2 emission of the CKs, particularly that produced in the interface 
H+ZH" of the knot, where a significant fraction of the H2 1-0 S(l) emission seems to 
be produced. Our results show that the production of molecular hydrogen in such a 
region may explain several characteristics of the observed emission, particularly the 
high excitation temperature of the H2 infrared lines. We find that the temperature 
derived from H2 observations even of a single knot, will depend very strongly on the 
observed transitions, with much higher temperatures derived from excited levels. We 
also proposed that the separation between the Ha and [Nil] peak emission observed 
in the images of CKs may be an effect of the distance of the knot from the star, since 
for knots farther from the central star the [N 11] line is produced closer to the border 
of the CK than Ha. 

Key words: Astrochemistry - circumstellar matter - ISM: molecule - infrared: ISM 
- planetary nebulae: individual: NGC 7293, Helix Nebula 



1 INTRODUCTION 

Planetary nebulae (PNe) are formed by the ejection of the 
outer layers of low and intermediate mass stars (1 to 8 Mq) 
in their final stages of evolution. The ejecta originate while 
the star is on the Asymptotic Giant Branch. Afterwards, the 
rapidly heating star drives a molecular dissociation front and 
an ionisation front, which eventually overrun the ejecta. The 
study of the composition and structure of PNe, in particu- 
lar, the molecular gas (its location, distribution, and phys- 
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ical conditions), can provide informat i on about the ea rlier 
evolutionary stages (|Davis et al.ll2003l : IVolk et al.ll2004l '). as 
well as on the phys ics of dissociation and ionisation fronts 
iHennev et al.ll2007l '). 

Infrared (IR) emission of H2 has been identified in mor e 
than 70 PNe jHora et al.ll 19991: ISter ling fc Dinerst"einll2008l '). 
primarily bipolar nebulae ( Kastner et al.lll996l '). The precise 
location of the observed H2 emission and the dominant ex- 
citation m echanism arc still under deb ate l|Hora et al.lll99sj : 
iLikkel et al,,2006, : Hcn nev et al.ll2007ll . 

Analysis of observations and models of the H2 IR line 
emission of PNe indicate that the emission is produced 
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both within the neutral envelope and within the ionised re- 
gion llAleman fc Gruenwald |2004| . and references therein). 
lAleman fc Gruenwaldl (2011) showed that, in cases of high 
temperature central stars, a significant part of the H2 emis- 
sion can be produced by the diffuse gas in the H+/H° tran- 
sition region of PNe. 

The Helix nebula (NGC 7293) is an evolve d PN excited 
Wa T* = 120000 K, = 100 L0 central star (|Henrv et all 
Il999l ). High-resolution images of this PN have shown that 
the H2 emission arises from its large population of glob- 
ules, the so-calle d cometary knots (CK s), embedded in the 
ionised gas (e.g. iMatsuura et al.ll2009( ). Many authors ar- 
gue that CKs are largely responsibl e for the H2 em ission 
produced within the ion i sed region (iBeckwith et al.l 
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Schildlll995l : ISpeck et aLll2002l : IMatsuura et al.ll2009t ) 

CKs are structures that resemble comets, particularly 
in images taken in Ha, [N 11] A6583, and H2 1-0 S(l) lines. 
Knot s and filamentary c ondensations are also seen in other 
PNe (lO'Dell et al.|[20o3 ) in the transition region between 
the ionised and neutral H. In the Helix nebula, the bright 
cusp points towards the central star and the tail in the op- 
posite direction, which can indicate that the excitation is 
connected with the central star. Analyses of the emission 
of these structures suggest that they are significantly denser 
than the gas around them. Whi le the typical dif fuse gas den- 
sity around 10^-10" cm"^ (|Osterbrock fc Ferland 200^ : 
lKwokll2000t), ill the CKs the density is arou r id lO'^-lO'^ cm~^ 
jHuggins et allll992l : iMeaburn et al.|[T998l : IMatsuura et al] 
|2007^ . The origin of the CKs is still uncertain. They may be 
formed within the original wind of the progenitor star or by 
instabilities and fragmentati on of the swept-up shell during 
the onset of the PN phase. iHuggins fc FranM (|2006l ) show 
that the latter is better supported by th e observed proper- 
ties of the knots. IMa tsuura et al.l (|2009l ) suggest a possible 
origin within a spiral density wave in the wind caused by an 
orbiting companion. 

There is no evi dence that the H2 emission in CKs is 
produced by s hocks ( Huggi ns et"aIll2002l : [O'Dell et al.ll2005l : 
IMatsuura et al. 200 7, 2008t) . On the other hand, as men- 
tioned above, the emission of th e CKs seems to be linked 
to the central star radiation field (|0'Dell et al.ll2005l . [20071 ). 
The H2 emission is more intense in a thin layer in the surface 
of the CKs towards the central star. However, analysis based 
on traditional models of photodissociation regions (PDRs) 
were unable to reproduce the high excitation temperature 
of the H2 emissi on (~900 K) estimated from the obser- 
vations (H uggins et all |2002| : lO'DeU et al.|[2007l ). Recently, 
iHennevet al.l (j2007 ) showed that advection can cause the 
ionisation and dissociation front to merge, leading to en- 
hanced heating of the molecular gas and reproducing well 
the excitation temperatures of H2. It is evident from the 
work of Henney et al. that the physical conditions in the in- 
terface between the ionized diffuse gas and the CK are of key 
importance for the H2 IR emission. The high excitation tem- 
perature of the H2 emission could be naturally explained if 
part of the H2 emission is produced in the H^/H" transition 
of the CK. 

In the present paper, the ionised and partially ionised 
regions of CKs are modelled with the one-dimensional pho- 
toionisation code Aangaba. Our aim is to study in more 
detail the H2 emission in the H+/H° interface of the CK. A 



grid of models is obtained to study how the H2 IR emission 
depends on different interface density profiles, CK radius, 
distance from the ionizing source, and dust-to-gas ratio. We 
compare the H2 IR emission with the atomic emission. Our 
calculations also provide good estimates of the H2 formation 
and destruction rates inside the CKs, which indicate how 
long the molecule may survive inside the CKs and whether 
it forms in-situ or may have survived the earlier evolutionary 
phase. 

We use our models to study the CKs in the Helix neb- 
ula (NGC 7293). The Helix is one of the nearest PNe (219 
pc), and high re solution i mages resolved the s tructure of 
the CKs (lO'De U ct al. 2004: Mcixncr e t al.ll2005l : lHora et all 
I2OO6I : IMatsuura et al.. 2007, . 2008. . . 2009 ). The detailed obser- 
vations allow us to test the model. Our models are described 
in Section 2. The results are discussed in Section 3 and con- 
clusions are summarized in Section 4. 



2 MODELS 

In Section 2.1, a brief description of the photoionisation code 
Aangaba is given. The parameters assumed for the Helix 
nebula model are discussed in Section 2.2. A description of 
how the CKs were simulated is given in Section 2.3. 



2.1 General Description of the Code 

The Aangaba code simulates the physical conditions in- 
side a photoionised nebula, given the ionising spectrum, as 
well as the characteristics a nd distribution of gas and dust 
(jCruenwal d fc Viegasl ll9"93 ). Since the gas and dust temper- 
atures and densities depend on each other, as well as on the 
position inside the nebula, the code makes iterative calcu- 
lations across the nebula. The code begins the calculations 
at the inner border of the nebula and continues the calcula- 
tion in steps in the outward direction. In each location, the 
physical conditions, as for example the ionic and molecu- 
lar density, the electronic density and temperature, the dust 
temperature, and continuum and line emission, are calcu- 
lated. The calculation ends when the code reaches the limit 
chosen by the user (for example, a given ionisation degree, 
gas temperature, optical depth, etc.). Geometrical dilution 
and extinction of the radiation by gas and dust are taken 
into account. The transfer of the primary and diffuse radia- 
tion fields is treated in the outward-only approximatiorQ. 



^ The outward-only is an usual approximation to treat 
the diffuse radiation in one-dimensional photoionization 
codes (see, for example, Pequignot et al. 2001; Stasinsl^ 
2009). It was introduced by Tarter & Salpeter (196^ 
In this approximation, in the words of Pequignot et alj 
(2001), 'all the diffuse photons are assumed to be emitted 
isotropically in the outer half-space'. The diffuse radia- 
tion photons emitted in the inward direction are assumed 
to travel to the point symmetrically opposite with re- 
spect to the central star in the other side of the nebula 
through the optically thin gas without being absorbed. 
This allows the code to proceed with the calculation in 
the outward direction only, that is, it does not need to 
recalculate previous steps to account for the diffuse ra- 
diation emitted more externally. 
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Twelve elements and their ions are included in the code: 
H, He, C, N, O, Mg, Ne, Si, S, Ar, CI, and Fe. The species 
H2, H^, H|^, and H~ are also includ ed in the code (H 3 
is not included because it is unstable; iBordas et a"l]|l99(J ). 
The density of each ion of each species is calculated with 
the assumption of chemical and ionisation equilibrium. The 
processes of ionization, recombination and charge exchange 
are taken into account in the equilibrium equations of the 
atomic species. For the H bearing species, forty reactions 
of formation and dest ruction are included in the chemical 
equilibrium equations jAleman fc Gruenwaldl [20041 '). 



The population of the H2 rovibrational levels of the 
three lowest electronic bound states is calculated by assum- 
ing statistical equilibrium, i.e., the total population rate of a 
level is equal to its total depopulation rate. For the electronic 
ground level, several excitation and de-excitation mecha- 
nisms (radiative and coUisional) , as well as the possibility 
that H2 is produced or destroyed in any given level, are in- 
cluded. For upper electronic levels, only radiative electronic 
transitions between each upper state and the ground state 
are included, since this must be the dominant mechanism. 
The population of the H2 rovibrational levels of the elec- 
tronic ground state by radiative mechanisms occurs through 
two main routes: electric quadrupole transitions between the 
rovibrational levels, involving IR photons, or electric dipole 
transitions to upper electronic states with subsequent decay 
to the ground state, involving UV photons (UV pumping). 
Collisions may also change the H2 energy level. We included 
collisions of H2 with the main components of the gas, i.e., 
H, H"'", He, H2, and electrons. In this work we are inter- 
ested in the lines produced by the rovibrational transitions 
of H2, whose wavelengths are in the 0.28 /xm to 6.2 mm 
range of the electromagnetic spectrum. More details about 
the calculation of the H2 lev el population can be found in 
lAleman fc Gruenwaldl l|201lh . 



The gas temperature is calculated assuming thermal 
equilibrium, that is, the total input of energy in the gas 
per unit time and volume is balanced by the total loss of 
energy per unit time and volume. Several mechanisms of 
gain and loss of energy by the gas due to atomic species, 
dust, and H2 are taken into account. The gas heating mech- 
anisms are photoionization of atoms, atomic ions, and H2 
by the primary and diffuse radiation; H2 photodissocia- 
tion (direct and two steps); H2 formation on grain sur- 
faces, by associative detachment, and by charge exchange 
with H; H2 coUisional de-excitation; and photoelectric ef- 
fect on dust surfaces, while the cooling mechanisms are 
emission of coUisionaly excited lines; radiative and dielec- 
tronic atomic recombination; thermal coUisional atomic ion- 
ization; free- free emission; coUisional excitation of H2; de- 
struction of H2 by charge exchange with H"*"; H2 coUisional 
dissociation; and co llision of gas-phase partic l es with dust 
grains. Models by 'Tielens HollenbachI (|l985l ') and 
[Spaans & Meiierink (2005) indicate other molecules 
become important contributors to the thermal bal- 
ance if Av > 3.5. Since in our models < 2.7, we do 
not include molecules other than H2 in the thermal 
equilibrium. 



Table 1. Input Parameters for the Helix Model. 





value 




120000 K 




100 Lq 


no 


50 cm-3 


Dust material 


amorphous carbon 


Grain radius 


0.1 /jm 


Distance 


219 pc 


Element 


Abundances (relative to H, by number) 


He 


1.20x10-1 





4.60x10-* 


C 


4.00x10-* 


N 


2.48x10-* 


Ne 


1.52x10-* 


S 


1.48x10-^ 


Ar 


3.10x10-^ 


Mg 


a.soxio-'^ 


Si 


3.50x10-'^ 


CI 


3.20x10-'' 


Fe 


4.70x10-'' 



References are given in the text. 



2.2 The Helix Model 



The parameters for the Helix nebula model are given in Ta- 
ble 1. We assume that the central star radiate s as a black- 
body with 7; = 120000 K and = 100 Lq (|Henrv et all 
I1999I : lO'DeU et ai]l2007l'). The dens ity of the diffuse gas is 
no = 50 cm-^ (|Meixner et~al]|2005h . 

Elemental abundances were determined bv lHenrv et al.l 

(|l999l ) for He, O, C, N, Ne, S, and Ar, from Helbc line 
emission observations and ICFs obtain ed from phq t onisa - 
tion models. The C/O ratio obtained bv lHenrv et al.l (|l999l ) 
is 0.87. The detection of C-bearing molecules, such as H2CO, 
C-C3H2, and C2H, has been used to argue th at the Helix may 
be in fact C-rich (|Tenenbaum et ah |2009|) . We ran some 
models with C-rich abundances, but we found it does not 
affect our results for H2 significantly, and will assume the 
iHenrv et all (| 19991 ) abundances. For the remaining elements 
taken into acco unt in the code, Mg, S i, CI, a nd Fe, we adopt 
the values f rom 'Stasihsk a fc Tvlendal (|l986l ). 

Spe ck et al. (2002) model the IR emission of Helix and 
found the dust-to-ga s ratio Md/Mg = 10-^, whi ch is an av- 
erage value for PNe (|Stasihska fc Szczer bal ll999l ). Our mod- 
els assume by default amorphous carbon grains with O.l/im 
radius. We also did calcul ations with silicate dust, but we 
found (as also discussed in ; Aleman & Grucnwald 2004) that 
the choice of compound does not cause significant changes 
in the H2 density. 

We assume a d istance of 219 pc ([Harris et al ] I2OO7I : 
[Benedict erai]|2009D . obtained from measurements of the 
central star parallax. 



2.3 Simulating the Cometary Knots 

The CKs are simulated as an increase in the density profile 
of the Helix nebula model at a given distance. We construct 
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Table 2. Parameters of the CKs models. 



Model 


Rk (arcsec) 


A/J (arcsec) 


(cm 3) 


Mi/M 


Kl 


129 


0.0 


10^ 


10-3 


K2 


129 


0.2 


10^ 


10-3 


K3 


129 


0.2 


IQ^ 


10-3 


K4 


129 


0.2 


IQ^ 


10-3 


K5 


129 


0.01 


10^ 


10-3 


K6 


129 


0.5 


10^ 


10-3 


K7 


129 


0.2 


10^ 


10-2 


K8 


129 


0.2 


lO'' 


10-3 


K9 


210 


0.2 


10^ 


10-3 


KIO 


383 


0.2 


10^ 


10-3 


Kll 


450 


0.2 


10^ 


10-3 


K12 


501 


0.2 


10^ 


10-3 



Profile " ^ 



RESULTS 
Diffuse Gas in the Helix Model 



a grid of CK models with different core densities (nx), den- 
sity profiles, dust-to-gas ratios (Md/Mg), and distances from 
the central star ( Rk)- We obtain models with uk between 
10^ to lO*^ cm-3 jHuggins et al.|[l99i : iMeaburn et"al]ll998l : 



iMatsuura et al.l 12007 ). We assume that the density profile 
has an increase from n_D to uk over a given distance A_R. 
In the following discussions, we call this region the interface 
of the knot. The region where the density reached the max- 
imum value (tik) is called the core. We study four types 
of density profiles: step function (hereafter type 1), linear 
(t ype 2), (type 3), and e xponential (type 4). According 
to lO'Dell fc Handrod l|l996t ). the Helix CKs have M^/Mg 
between 7x10-* and 7x10-^. We obtain models for Md/Mg 
of 10-3 and IQ-'^. In each model, the dust-to-gas ratio and 
the chemical composition are assumed the same for CKs and 
diffuse gas. The parameters of the CK models discussed in 
the present work are listed in Table [2] 

The code starts the calculations in the inner border of 
the nebula and continues outward along the radial direction. 
At a given distance (Rk), an increase in density simulates 
the CK. For the present work, we calculate models with 
Rk = 130 to 500 arcsec (see discussion in Sect. 13.1] ). We stop 
the calculations where the gas temperature, which decreases 
with distance from the central star, reaches 100 K. We define 
ARt as the distance from the border of the knot to this 
point. 

An IDL routine was developed to simulate a three- 
dimensional CK, allowing the calculation of surface bright- 
ness of atomic and H2 lines by the integration of the emis- 
sivity along the line of sight inside the CK. In this routine, 
the knot is assumed to have cylindrical symmetry, with a 
semi-spherical head pointing in the direction of the central 
star. The CK is assumed to be seen edge-on. We assume 
that the one-dimensional profile of the emissivity, which is 
calculated with Aangaba, is the same for every line along 
the direction parallel to the symmetry axis, starting from 
the border of the knot. This is a simple approximation to 
simulate a comet-shaped knot, but it allows us to obtain a 
reasonable estimate of the H2 lines surface brightness. It is 
not our intention to reproduce the precise image of a CK. 



Tigure [T] shows results for the Helix model obtained with 
the parameters of Table [T] The gas temperature profile is 
given in the top panel, while the radial profiles of H'^, H"^, 
H-, HJ, and HJ densities are plotted in the middle panel. 
Emissivities of some lines are shown in the bottom panel. 
The model in the figure applies to the diffuse gas only: it 
does not include CKs. The gas density and the dust-to-gas 
ratio are assumed to be uniform and equal to, respectively, 
ud = 50 cm-3 Md/Mg = 10-3. 

In the Helix, CKs are detected at distances from the 
central star (pr ojected on the sky plan e) between 132 arcsec 
.and 384 arcsec l|Matsuura et al.ll2009l ). This interval is indi- 
cated by the gray band in Fig. [T] The three-dimensional 
morphology must be taken into account to calculate the 
real distance. The m orph ology of the Helix w as studied by 
I and lO'Dell et all l|2004h . According to 
^), the helix structure (where the CKs 
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are) has a toroidal shape, contained within an angle of about 
20°on both sides of the equatorial plane of the nebula. The 
symmetry axi s is inclined a t abo ut 37°with respect to the 
line of sight. lO'DeU et all l|2004l ) found a shght different 
value to the inclination angle, 28° with respect to the line 
of sight. Taking this geometry into account, real distances 
must be greater than the projected distances. For example, 
a CK in a plane inclined by 37° from the line of sight could 
be up to 25 per cent farther from the central star than the 
projected distance. Even with this deprojection, the CKs are 
not much farther than 500 arcsec from the central star, that 
is, mostly inside the region where hydrogen is ionised. Our 
model also supports that the CKs coexist with the ionised 
species He"*", N^,N^"'", O^, and 0^"*" in the diffuse gas. Our 
model has parame ters similar to those of model 135/315 of 
iHenrv et al.l (Il999l ) and provides similar results for the ion- 

isation structure. 

According to IMatsuura et aL I l|2009l ) andlM 

eixner et al.l 

the H2 emission of the Helix is associated with the 
CKs, although the instrumental resolution can not exclude 
a contribution from the diffuse gas in the outer ring. As 
can be seen in Fig. [1] the calculated H2 1-0 S(l) emissiv- 
ity from the diffuse gas is very low. Assuming a spheri- 
cal gas distribution, the maximum H2 1-0 S(l) intensity is 
10"** erg cm-2 s-'^ sr-^, and the total 1-0 S(l) flux com- 
ing from inside the ionised region is 4x10-^3 erg cm^^ g-i^ 
In contrast, the observed 1-0 S(l) flux of the Helix nebula 
is about 2x10-^ erg cm-^ s-^ (obtai ned from the obser- 
vations published in ISpeck et al. I I2OO3I ). The flux obtained 
with our model without GCs shows that the diffuse ionised 
gas contributes less than 0.02 per cent of the observed emis- 
sion. Our calculation of the flux above does not take into 
account a possible increase in the average gas density in the 
rings or the shielding of the radiation by the CKs, which 
may increase the H2 emission. 



3.2 Physical and Chemical Properties of CKs 

Figures [21 131 and (31 show results for models Kl to K12 (see 
the model parameters in Table [2]| . The upper panels show 
the gas temperature and the total H nuclei density profile; 
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Figure 1. Diffuse gas in tlie Helix model: (a) gas temperature, 
(b) density of H", H+, H~, H2, H^, and H^; and (c) emissivity 
of some lines. CKs are not included in this model. The grey band 
represents the range of distances where the CKs are observed in 
the Helix; see discussion in the text. 



the middle panels show the density profiles of H", H^, H , 
H2, HJ, and Hj^; and the bottom panels show the emissivity 
profiles of H2 1-0 S(l), Ha, [Nil] A6583, [Olll] A5007, and 
[Ol] 63^111. 

The gas temperature decreases with depth from the lo- 
cal diffuse gas value at the border of the CK to the 100 K 
limit assumed for our calculations. In the interface, the gas 
is hot, with temperatures T > 2000 K. As previously men- 
tioned, the balance between the energy gain and loss by the 
gcis determine its temperature. Figure [5] shows the rate of 
energy gain and loss inside the CK model K4. The results 
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Figure 5. Contribution of the mechanisms involving atoms, dust, 
and H2 to the rates of energy gain and loss by the gas inside the 
CK model K4. 



are analogue for other models. In the interface, atomic pro- 
cesses dominate the heating and cooling of the gas, as in the 
diffuse ionised gas around the CK. The H2 molecule domi- 
nates the heating and cooling in the core of the CK, while 
atomic mechanisms contributes to about 10 per cent of the 
total rates of energy loss and gain. CoUisional excitation and 
de-excitation of H2 are the most important mechanisms of 
loss and gain, respectively. Photoionisation of II2 may also 
contribute to the energy gain (<10 per cent). Heating by 
the photoelectric effect on grains may be significant only if 
Md/Mg » 10"^. For example, this mechanism contributes 
less than one per cent to the total energy input in model 
K4 (Md/Mg = 10~^), and less than 5 per cent for model K7 
(Md/Mg = 10"^). The calculated grain temperature inside 
the CKs is approximately 20 K. 

The transition between ionised and neutral H occurs in 
the interface, around T ~ 8000 K. Below this temperature 
and down to 100 K, hydrogen is mostly neutral. The den- 
sity of molecular hydrogen increases inward within the CK, 
reaching a fractional abundance in excess of 30 per cent. In 
cases with a high dust-to-gas ratio (see Sect. l3.2~T1l . a fully 
molecular region is produced while the temperature is still 
above 100 K. For other dust-to-gas ratios, the fully molecular 
region must have lower temperatures. The relative densities 
of H~, HJ, and Hj are smaller than 10~^ in all the models. 

According to our models, inside a CK the most impor- 
tant H2 formation mechanisms are: 
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Figure 2. Results for models Kl to K4 (left to right). Top: gas temperature (dot-dashed) and total density (solid) profiles. Middle; 
density of H" (dashed), H+ (dot-dashed), H~ (long dash), H2 (solid), (three dot-dashed), and Hj{" (dotted). Bottom: cmissivities of 
the lines H2 1-0 S(l) (soUd), H« (dashed), [Nil] A6583 (dot-dashed), [Olll] A5007 (long dashed), and [Ol] G^i^^ul (dotted). Distance is 
taken from the border of the knot facing the star along the symmetry axis. 



(Fl) H -f H- ^ H2 e- 
(F2) H+ H ^ H+ -f- Ha 
(F3) H+ e- ^ Ha -f H 
(F4) 2H + grain ^ Ha + grain 

and the main Ha destruction mechanisms are: 

(Dl) Ha + hi/ -> H+ 6" 
(D2) Ha + hi/ -> 2H 

(D3) Ha + hi/ ^ Ha 2H 

(D4) Ha + H+ ^ H H+ 

(D5) Ha + H+ ^ H+ H 
(D6) Ha + H ^ 3H 

Both neutral and ionised species are involved in the Ha for- 
mation and destruction processes; UV photons are impor- 
tant for the destruction of the molecule. Figure [6] shows the 
rate of these processes as a function of the distance inside 
a CK. The associative detachment (Fl) and the charge ex- 
change (F2) reactions dominate the formation of Ha at the 
border of the CK facing the star. The associative detachment 
is less important towards high depths, while the formation 
on grain surfaces becomes more important and dominates 
the Ha formation in the CK core. The charge exchange and 
ion-molecule (F3) reactions contribute significantly to the 
molecular formation in the whole studied region. 



The list of important reactions of Ha formation and de- 
struction in CKs is si milar to the im portant reactions for the 
diffuse ionised gas (see I Aleman fc G rucnwald 2004). The dif- 
ferences are: (F3) and (D5) are not important in the diffuse 
gas, and the dissociation of Ha by electron collision is not 
important insid e the CKs. 

Models by iHennev et al.l |20o3) provide almost the 
same list of significant reactions for Ha, with the exception 
of destruction of Ha by the reaction Ha + e^ — >■ H^ + H, 
instead of (D6), probably because in their models the ad- 
vection brings the ionisation front closer to the molecular 
zone, increasing the density of electrons and decreasing the 
density of H'^ coexisting with Ha. 

Estimates of the timescale for Ha formation/destruction 
can be obtained with the formula 



dn{H2)/dt' ^> 

where At is the characteristic timescale, n(Ha) is the 
molecular hydrogen density, and An /At is the for- 
mation/destruction rate. Our equilibrium models give 
timescales inside a CK typically from one year at the border 
to 10'* years in the core of a CK. 
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Figure 3. Same as Fig. [2] for models K5 to K8 (left to right). Note the different density scale for model K8 in panel (a8). 



3.2.1 Effect of the Input Parameters 

As shown in Figs.[2]to|4l models with different input param- 
eters yield different ionisation structures and hence result 
in different line emissivity profiles. Figure [2] shows models 
with different CK density profiles (models Kl to K4). Us- 
ing the parameters of the model K4 as reference, in model 
K5 and K6 we change the interface thickness; in model K7, 
the dust-to-gas ratio; in model K8, the core density (Fig. [3]). 
Models K9 to K12 have different distances from the central 
star (Fig.[3)). Figure [7] summarises the effect of the input pa- 
rameters on the total thickness from the border of the knot 
to the position where T = 100 K (i.e., A7?t). The figure 
also shows the thickness of ionised and neutral zones of this 
region. Here the boundary between the ionised and neutral 
zones is assumed to be where H ionization degree is 0.5. 

The thickness of the CK ionised zone increases signif- 
icantly from type 1 to type 4 interface. The variation of 
the thickness of the neutral zone is less significant, but it 
also changes with the type of the interface. The model with 
type 1 interface has the smallest absolute thickness, while 
type 2 has the greatest value. As a result, A_Rt increases 
significantly from type 1 to type 4 interface. 

The density profile in the interface may result from the 
physical processes occurring and the nature of knots. For 
example, an function woul d be expected in the cas e of a 
spherical e quilibrium outfiow jBurkert fc 0'Dellll998l ). Ac- 
cording to iBurkert fc O'Delll (|l998l ') the Hq surface bright- 



ness profile is better fitted by an exponential function than 
by an r'^ function (resu lting from an r^ increase in density). 
Following the results of lBurkert fc O'Delll ()l998h . we assume 
an exponential increase hereafter. 

The thickness of the ionised zone and Ai?T increase 
for models with thicker interfaces (greater A_R), as can be 
seen by comparing models K5, K4, and K6 (Figs. [2] and 
[31 . The three models have the same parameters, except for 
Ai?, which takes, respectively, the values 0.01, 0.2, and 0.5. 
In models with Ai? = 0.01 the neutral region is about 10 
times larger than the ionised, while for AJ? = 0.2 they are 
about the same size. 

The rate of H2 formation on grain surface is propor- 
tional to the dust-to-gas ratio. The contribution of dust to 
the gas energy input also increases in models with higher 
Md/Mg, producing a more extended region with 100 < T < 
300 K. This allows hydrogen to become fully molecular in 
such a region. Compare models K4 (Fig. [2)) and K7 (Fig. 
[Sjl, which differ only by the value of dust-to-gas ratio (the 
dust-to-gas ratio for K7 is 10 times larger than for K4). 

Density and gas temperature profiles for model K8 is 
shown in Fig. [3] This model is similar to model K4, but with 
a higher core density tik. Comparing both models, it can be 
seen that both ionised and neutral regions reduce in size in 
the denser model. Several mechanisms can contribute to this 
difference, such as the increase in the radiation extinction by 
gas and dust, increase in the H2 formation on grain surface. 
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Figure 4. Same as Fig.Elfor models K9 to K12 (left to right). 



heating of the gas. It is important to notice that, since we 
keep the dust-to-gas ratio constant throughout the model, 
an increase in the gas density implies in an increase of the 
dust density. 

An important parameter to the ionisation structure of 
the CKs is the distance from the central star, since the ion- 
ising flux and spectrum may change significantly with the 
position in the nebula. CKs farther from the central star 
have smaller ionised zones. If the CK is beyond the Helix 
ionisation front there is no ionised region, and the knot is 
completely neutral (see model K12 in Fig. [4]). 



3.3 Warm H2 1-0 S(l) Emission 

As can be seen in Figs. [2] to [4] the emissivity of the H2 1- 
S(l) line in the CKs is high in a warm region, where the 
temperature ranges between 300 and 7000 K. The peak in 
the 1-0 S(l) emissivity in the studied region occurs where the 
density is around 40 per cent of the core density (with the 
exception of interface type 1 models, which have no interme- 
diate density). This is also true for other rovibrational lines. 
This component of the H2 emission can explain the excita- 
tion t empe ratures around 9 00-18 00 K found by ICox et al] 
l|l998h and lMatsuura et all |20o3). The spatial width of the 
peak in the H2 1-0 S(l) emissivity is similar for all models, 
but the peak emissivity decreases with the distance from the 
central star. 



Inside the CKs, the H2 rovibrational levels of the ground 
electronic state are mainly excited and de-excited by coUi- 
sional processes, particularly for lower J levels. Formation 
pumping may contribute to high J levels population (J > 
10) , in the region where H is neutral. In such a region, rovi- 
brational radiative de-excitation is an important mechanism 
for the population of levels with J > 15 or for lower J levels 
with w ~ 5. UV Pumping is important for the population of 
levels with w > 4 and lower J. 

Figure|H]shows the H2 excitation diagram for model K9. 
The effective column densities of the rovibrational levels v — 
0, 1, 2 and 3 were calculated using the intensities determined 
in our model: these are indicated by the filled symbols. The 
lines in Fig. [8] represent Boltzmann distributions for three 
different temperatures. Different regions within the CK con- 
tribute to different parts of this plot. The highest energy lev- 
els are populated by collisions in the hottest regions of the 
CK. The lower energy levels are populated throughout the 
knot, at a range of temperatures. As a result, the integrated 
spectrum yields a much higher temperature for transitions 
involving a high upper energy level. The plot clearly shows 
these different excitation components of the II2 emission. 
The lines of the band 1-0 and 2-1 are well represented by 
an excitation temperature of approximately 2000 K. The 
0-0 band shows a distribution of lower temperatures, ap- 
proximated in the plot as two components, with excitation 
temperatures of 900 and 200 K. 
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Figure 6. Rates of the important processes of H2 formation and 
destruction inside the CK model K4. See text for the reaction 
codes. 



Effective col umn densities obtained from peak intensi- 
ties measured bv lMatsuura et al] l|2007t ) (their region a) are 
shown in Fig. [8] (open symbols). The agreement between the 
excitation temperatures of t he m odel and the observations 
is evident. iMatsuura et al.l (|2007l ) found a temperature of 
1800 K from the 1-0 and 2-1 bands in a knot relatively close 
to the star, well represented by the model for these bands. 
(Note that the observed knot is 150 arcsec from the star, 
while the model is for a (deprojected) distance of 210 arc- 



For the 0-0 band. lCox eFal] (|l998l ') found from ISO ob- 
servations that the observed column densities obtained from 
lines S(2) to S(7) are well represented by Boltzmann distri- 
bution at a temperature around 900 K. This agrees with the 
excitation temperature of 1000 K predicted by the current 
model. 

IMatsuura et al.l (|2007l ) suggested that the different exci- 
tation temperatures they found (inferred from ro vibrational 
lines) and the value obtained bv lCox et al.l (|l998t ) (from pure 
rotational li nes) is evidence for temperature variation within 
the nebula. iHora et alj (|200 6f) found some indication for a 
gradient in the excitation in the Helix CKs with the distance 
from the central star, but in our models, a significant gradi- 
ent is expected only at large distances, close to the ionisation 
boundary of the PN . Instead, t he dif ferent excitation tem- 
peratures found by I Cox et all l| 19981 ) and IMatsuura et al.l 




200 300 400 500 
Rk (arcsec 

Figure 7. Effect of the various input parameters on the H ioni- 
sation structure. In each plot, the total thickness from the border 
of the knot to the position where T = 100 K (ARt) is repre- 
sented by the solid line. The widths of the ionised and neutral 
zones are represented by the dot-dashed and the dashed curves, 
respectively. 



|2003) can be explained by the effect of the excitation within 
a single CK. 

Using the ratios between lines of the S and Q branches 
of band 1-0 lines departing from the same level and assumin g 
an extinction law of the form A~ (see iDavis et al.ll2003l ). 
we estimate an extinction of about 10-20 per cent at 2/^m. 
The model assumes that the knot is fully transparent, while 
in reality a fraction of the emission from the backside of the 
knot is extincted. However, this will not have a major effect 
on our results. 

For comparison with H2, the emissivity of a few in- 
tense atomic lines are also shown in Figs. [2]to |3| The [Nil] 
A6583 emission is produced in the interface of the CKs. Ha 
is also produced in the interface region, but there is a less in- 
tense extended emission region produced by the high energy 
photons (which ionise the H atoms; the following recombi- 
nation produces the line). The emissivity of the [Nil] line 
loses importance for models farther from the central star. 
For such models, the [Nil] line is produced only near the 
border of the CK and its emissivity profile becomes very 
different from that of Ha. Such behaviour may explain the 
separation between the Ha and [Nil] peaks observed by 
lO'Dell et ai] (|2000l ).The decrease with distance of the [Nil] 
emissivity is stronger than for the 1-0 S(l) line, which can 
explain why all CKs with [Nil] emission have a 1-0 S(l) 
counterpart, wh i lst th e opposite is not true, as observed by 
IMatsuura et al.l l|2009l '). 

Inside the CKs, the [O lll] A5007 emission is only pro- 
duced in the interface and in models closer to the central 
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Figure 8. H2 excitation diagram. The effective column density 
was calculated from the surface brightness at the peak intensity 
of 1-0 S(l). Open symbols represent observations and filled sym- 
bols represent models. Lines are Boltzmann distributions for the 
temperature indicated near each curve. 



star. In such cases, the emissivity inside the CK increases 
only up to one order of magnitude from the diffuse gas value, 
a small value when compared to other lines, such as Ha or 
[N 11] A6583. Images of the Helix show [O ill] emissio n in the 
diffus e gas and, in some cases, in tlic tip of CKs ( O' Dell et al.l 
bOOTh . According to IWalsh fc MeaburnI l|l993f ). 40 per cent 
of the CKs does not show any emission in [Olll] images. 

The predicted emissivity of the IR fine structure line 
[O l] 63^m inside the CKs is very high as can be seen in 
Figs. [2j[3l and|4] Reay et al (1988) found a correlation be- 
tween the H2 1-0 S(l) and the optical [OI] A6300 in a sam- 
ple with several PNe, including the Helix. The optical line 
is produced only in the interface, while the profile of the IR 
[01] line is more similar to the one of 1-0 S(l). 

As discussed before, if the CK is beyond the Helix ioni- 
sation front there is no ionised region, and the knot is com- 
pletely neutral (see models Kll and K12 in Fig. |4}. In this 
case, the emissivity of both [N 11] and [O l] lines mentioned 
above is very low in the CKs and would not be detected. 




100 



200 300 
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Figure 9. H2 1-0 S(l) surface brightness of a cometary knot 
as a function of the distance to the central star. Sets of solid, 
dashed, and dot-dashed curves represent models with Ai? = 0.5, 
0.2, and 0.01 arcsec, respectively. Different curves within each set 
represents CK radius of 0.5, 1.0, and 2.0 arcsec, with the sur- 
face brightness increasing for larger CK radius. Dots represent 
measured values. The uncertainty in the distance is estimated as- 
suming that the Helix symmetry axis is inclined 37 ° with respect 
to the line of sight 



Table 3. H2 1-0 S(l) Surface Brightness Measurements. 



RA (E2000) <5 (E2000) Distance Peak Surface Brightness 



(h m s) 


(°"') 


(arcsec) 


(erg cm ^ 


s-l 


sterd -*-) 


22:29:45.987 


-20:49:05.98 


125 


2.0 


xlO 


-5 


22:29:42.852 


-20:49:01.39 


95 


3.7 


xlO 


-5 


22:29:42.873 


-20:47:43.35 
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1.7 


xlO 


-5 


22:29:47.370 


-20:47:37.92 


199 


1.2 


XlO 


-5 


22:29:40.616 


-20:47:52.70 


144 


9.6 


XlO 


-6 


22:29:41.285 


-20:47:14.67 


183 


9.2 


XlO 


-6 


22:29:37.829 


-20:48:01.66 


132 


2.8 


XlO 


-5 


22:29:42.820 


-20:48:27.77 


122 


1.3 


XlO 


-5 


22:29:39.823 


-20:47:53.55 


141 


1.9 


XlO 


-5 



3.4 H2 1-0 S(l) Surface Brightness 

Measurements of the H2 1-0 S(l) line surface brightness of 
some representative CKs in the Helix nebula are plotted in 
Fig. |9l We identified 10 isolated CKs detected both in Ho 
and H2 images. We m easured 2.12 ^m H2 in tensities from 
the images obtained bv lMatsuura et all (|2009l ). To calibrate 
the intensities, we used five stars within the observed field 
to measure the zero-point. We assume that the 2MASS K'- 
magnitude of these stars are the same as the magnitudes 
in H2 filter. We apply a 25-pixel radius for aperture pho- 
tometry and take the 35-50 pixel ring for the background 
measurements. The pixel scale is 0.117 arcsec. Table [3] lists 
the knots positions and the measurements. 

Curves in Figure |9] show the calculated H2 1-0 S(l) sur- 
face brightness of CK models as a function of the CK dis- 
tance to the central star. The calculated surface brightness 
was averaged over the same aperture as the measurements 



to allow direct comparison. Different curves represent dif- 
ferent AJ? and CK radius. The surface brightness increases 
with a decrease of Aii, an increase of the CK radius. CKs 
closer to the central star tend to have higher 1-0 S(l) surface 
brightness. The decrease in the H2 1-0 S(l) surface bright- 
ness with dis t ance from the central star was also noticed by 
lO'Dell et ai] l|2007h . If the CK is beyond the Helix ionisa- 
tion front, the 1-0 S(l) surface brightness drops dramati- 
cally, since there would be no enough radiation or temper- 
ature to excited significantly the upper vibrational levels of 
the molecule. Models can reproduce the magnitude of ob- 
served surface brightness and its decrease with distance to 
the central star. 

The 1-0 S(l) peak brightness is slightly higher for mod- 
els with higher dust-to-gas ratio. The increase caused by 
changing the dust-to-gas ratio from 10~^ to 10~^ is about 
20 per cent. The 1-0 S(l) peak brightness also increases with 
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increasing uk- The difference between models with uk = 
10* and 10^ cm~"^ is up to 40 per cent (20 per cent in the 
models farth er from the ce ntral star). 

.Matsuura et al. I l|2009l ) show that, for the same number 
of CKs, the outer ring shows twice the surface brightness 
as the inner ring. According to them, such difference may 
be caused by the diffuse gas or by the difficulty in isolat- 
ing CKs in the images. Since we do not expect a significant 
contribution from the diffuse gas and the surface brightness 
decreases with the distance to central star, the second ex- 
planation seems more likely. 



4 CONCLUSION 

As discussed in the Introduction some features of the CKs 
emission indicate that they are significantly affected by the 
radiation field produced by the Helix central star. In this 
paper we use the photoionisation code Aangaba to as- 
sess the contribution that photoionisation may have on the 
molecular hydrogen emission, supporting that photoionisa- 
tion plays a major role in such o bjects, which i s in a gree- 
ment with previous results from lO'Dell et al.l (|2007l ) and 
iHennev et ahl l|2007f ). 

Our Helix nebula model with no CKs indicates that 
emission of H2 from ionised diffuse gas is not significant, 
agreeing with the observations that the H2 emission of the 
Helix is associated with the CKs dMatsuura et al.l l2009l : 
iMeixner et "ai]|2005l ') 

We study CK models with different sets of input pa- 
rameters, chosen to cover the range of values inferred from 
Helix observations. The important mechanisms of gas heat- 
ing and cooling, formation and destruction of H2, excitation 
and de-excitation of H2 were determined. 

The emissivity of the H2 1-0 S(l) line in the CKs is 
important in a region with temperatures between 300 and 
7000 K. Such a warm component of the H2 emission may 
explain t he excitation tem per atures of around 900-1 800 K 
found bv lCox et all ll998') and ' Matsuura et alj (|2007l '). The 
contribution of this warm region to the emission of rovibra- 
tional lines is very important. For pure rotational lines of the 
II = level the contribution of the colder regions should be 
significant. The excitation diagrams obtained from our mod- 
els seems to agree very well with the observations. Compari- 
son with measurements indicates that models can reproduce 
the magnitude of observed surface brightness. 

An important parameter to the ionisation structure of 
the CKs is the distance from the central star, since the ion- 
ising spectrum may change significantly with the position in 
the nebula. CKs farther from the central star have smaller 
ionised zones. If the CK is beyond the Helix ionisation front 
there is no ionised region, the knot is completely neutral 
and there would be no enough radiation or temperature 
to excited significantly the upper vibrational levels of the 
molecule. The 1-0 S(l) intensity would be very low. 

Our results also indicates tha t the separation bet ween 
the Ha and [N 11] peaks observed bv lO'Dell et al.l (|2000t ) may 
be an effect of the distance of the knot from the star, since for 
knots farther from the central star the [N 11] line is produced 
closer to the border of the CK than Ha. 

The decrease with distance of the [Nil] emissivity is 
stronger than for 1-0 S(l) line, which can explain why all 



CKs with [N 11] emission have a 1-0 S (l) counterpart, but th e 
opposite is not true, as observed by Matsuura et al.l |200^. 

As pointed out by ' Burkert fc O'Delll IiQqI T the inter- 
face between the diffuse gas and the CK core may provide 
important clues about the mechanisms that shape and sus- 
tain the PNe CKs. Our models show that there are sig- 
nificant differences in the models results depending on the 
assumed density profiles of this region. Images that could 
resolve this region in great detail are then essential. 

We find that the temperature derived from H2 observa- 
tions even of a single knot, will depend very strongly on 
the observed transitions, with much higher temperatures 
derived from excited levels. This is caused by the large 
range of temperatures present, and the presence of signif- 
icant amounts of molecular hydrogen within the mini-PDR. 
This explains the puzzling temperature differences found by 
previous observers. 

We also find that H2 required to fit the observations, is 
consistent with the abundances calculated through equilib- 
rium chemistry, with a variety of formation and dissociation 
reactions. There is no need to presume that the molecular 
hydrogen in the knots predate the ionisation of the neb- 
ula, assuming that there is enough time to reach equilib- 
rium. This does not resolve the issue of when and how the 
knots formed, but we cannot exclude models where the knots 
formed after the onset of ionisation, nor models where the 
knots formed already in the AGB wind. 



ACKNOWLEDGMENTS 

LA. acknowledged the financial support from CNPq Brazil 
PDE fellowship number 201950/2008-1 and STFC rolling 
grant ST/F003196. LA. is grateful to the staff and students 
of the Astrophysics Group of University College of London 
for the hospitality during her visit in May 2010. We thank 
Paul Woods for discussions about the chemistry of C and O 
rich planetary nebulae, and the anonymous referee for the 
valuable suggestions to improve this paper. 



REFERENCES 

Aleman I., Gruenwald R., 2004, ApJ, 607, 865 
Aleman I., Gruenwald R., 2011, A&A, 528, A74H- 
Beckwith S., Gatley I., Persson S. E., 1978, ApJ, 219, L33 
Benedict G. F., McArthur B. E., Napiwotzki R., Harrison 

T. E., Harris H. C, Nelan E., Bond H. E., Patterson R. J., 

CiarduUo R., 2009, AJ, 138, 1969 
Bordas C, Cosby P. C, Helm H., 1990, J. Chem. Phys., 

93, 6303 

Burkert A., O'DeU C. R., 1998, ApJ, 503, 792 

Cox P., Boulanger F., Huggins P. J., Tielens A. G. G. . M., 

Forveille T., Bachiller R., Cesarsky D., Jones A. P., Young 

K., Roelfsema P. R., Cernicharo J., 1998, ApJ, 495, L23-I- 
Davis C. J., Smith M. D., Stern L., Kerr T. H., Chiar J. E., 

2003, MNRAS, 344, 262 
Gruenwald R. B., Viegas S. M., 1992, ApJS, 78, 153 
Gussie G., Pritchet C, 1988, JRASC, 82, 69 
Harris H. C, Dahn C. C, Canzian B., Guetter H. H., 

Leggett S. K., Levine S. E., Luginbuhl C. B., Monet 



12 /. Aleman et al. 



A. K. B., Monet D. G., Pier J. R., Stone R. C, Tille- 

man T., Vrba F. J., Walker R. L., 2007, AJ, 133, 631 
Henney W. J., Williams R. J. R., Ferland G. J., Shaw G., 

O'Dell C. R., 2007, ApJ, 671, L137 
Henry R. B. C, Kwitter K. B., Dufour R. J., 1999, ApJ, 

517, 782 

Hora J. L., Latter W. B., Deutsch L. K., 1999, AJS, 124, 
195 

Hora J. L., Latter W. B., Smith H. A., Marengo M., 2006, 
ApJ, 652, 426 

Huggins P. J., Bachiller R., Cox R, Forveille T., 1992, ApJ, 
401, L43 

Huggins P. J., Forveille T., Bachiller R., Cox P., Ageorges 

N., Walsh J. R., 2002, ApJ, 573, L55 
Huggins P. J., Frank A., 2006, in M. J. Barlow & 

R. H. Mendez ed.. Planetary Nebulae in our Galaocy and 

Beyond Vol. 234 of lAU Symposium, The formation of 

globules in planetary nebulae, pp 271-276 
Kastner J. H., Weintraub D. A., Gatley I., Merrill K. M., 

Probst R. G., 1996, ApJ, 462, 777 
Kwok S., 2000, The Origin and Evolution of Planetary 

Nebulae. Cambridge ; New York : Cambridge University 

Press, 2000. (Cambridge astrophysics series ; 33) 
Likkel L., Dinerstein H. L., Lester D. F., Kindt A., Bartig 

K., 2006, AJ, 131, 1515 
Matsuura M., Speck A., Smith M., Zijlstra A., Lowe K., 

Viti S., Redman M., Wareing C, Lagadec E., 2008, The 

Messenger, 132, 37 
Matsuura M., Speck A. K., McHunu B. M., Tanaka I., 

Wright N. J., Smith M. D., Zijlstra A. A., Viti S., Wesson 

R., 2009, ApJ, 700, 1067 
Matsuura M., Speck A. K., Smith M. D., Zijlstra A. A., Viti 

S., Lowe K. T. E., Redman M., Wareing C. J., Lagadec 

E., 2007, MNRAS, 382, 1447 
Meaburn J., Clayton C. A., Bryce M., Walsh J. R., Hol- 

loway A. J., Steffen W., 1998, MNRAS, 294, 201 
Meixner M., McCuUough P., Hartrnan J., Son M., Speck 

A., 2005, AJ, 130, 1784 
O'Dell C. R., Balick B., Hajian A. R., Henney W. J., Burk- 

ert A., 2002, AJ, 123, 3329 
O'Dell C. R., Handron K. D., 1996, AJ, 111, 1630 
O'Dell C. R., Henney W. J., Burkert A., 2000, AJ, 119, 

2910 

O'Dell C. R., Henney W. J., Ferland G. J., 2005, AJ, 130, 
172 

O'Dell C. R., Henney W. J., Ferland G. J., 2007, AJ, 133, 
2343 

O'Dell C. R., McCuUough P. R., Meixner M., 2004, AJ, 

128, 2339 

Osterbrock D. E., Ferland G. J., 2006, Astrophysics of 
gaseous nebulae and active galactic imclei. 2nd. ed. by 
D.E. Osterbrock and G.J. Ferland. Sausalito, CA: Uni- 
versity Science Books, 2006 

Pequignot D., Ferland G., Netzer H., Kallman T., Ballan- 
tyne D. R., Dumont A., Ercolano B., Harrington P., Krae- 
mer S., Morisset C, Nayakshin S., Rubin R. H., Suther- 
land R., 2001, in G. Ferland & D. W. Savin ed.. Spec- 
troscopic Challenges of Photoionized Plasmas Vol. 247 
of Astronomical Society of the Pacific Conference Series, 
Photoionization Model Nebulae, pp 533 — h 

Reay N. K., Walton N. A., Atherton P. D., 1988, MNRAS, 
232, 615 



Schild H., 1995, A&A, 297, 246 

Spaans M., Meijerink R., 2005, ApSS, 295, 239 

Speck A. K., Meixner M., Fong D., McCuUough P. R., 

Moser D. E., Ueta T., 2002, AJ, 123, 346 
Speck A. K., Meixner M., Jacoby G. H., Knezek P. M., 

2003, PASP, 115, 170 
Stasinska G., 2009, What can emission lines tell us?, pp 1- 

+ 

Stasinska G., Szczerba R., 1999, A&A, 352, 297 
Stasinska G., Tylenda R., 1986, A&A, 155, 137 
Sterling N. C, Dinerstein H. L., 2008, AJS, 174, 158 
Tarter C. B., Salpeter E. E., 1969, ApJ, 156, 953 
Tenenbaum E. D., Milam S. N., Woolf N. J., Ziurys L. M., 

2009, ApJ, 704, L108 
Tielens A. G. G. M., 1993, in Weinberger R., Acker A., eds. 
Planetary Nebulae Vol. 155 of lAU Symposium, Photodis- 
sociation Regions and Planetary Nebulae, pp 155 — h 
Tielens A. G. G. M., HoUenbach D., 1985, ApJ, 291, 722 
Volk K., Hrivnak B. J., Kwok S., 2004, ApJ, 616, 1181 
Walsh J. R., Meaburn J., 1993, The Messenger, 73, 35 

This paper has been typeset from a TisX/ X^T^jK file prepared 
by the author. 



